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The aim of this paper is to evaluate different types of kinetic models and to propose a set of kinetic 
parameters for simulating the heat treatment of woody biomass under temperatures representative 
of torrefaction conditions. The mass loss of three wood types, locust, spruce and eucalyptus, as well 
as their main components, cellulose, xylan, glucomannan and lignin, has been determined in a TGA at 
different isothermal temperatures between 210 and 290 °C for residence times of up to ten hours. These 
data were used to evaluate the kinetic parameters corresponding to various multi-step series reaction 
mechanisms proposed in the literature. A computational solution of the multi-step formulation was 
derived, which allows the exact conditions undergone by the sample to be used as input parameters. The 
inverse procedure takes into account simultaneously the different temperature levels for each species. 
A comprehensive table summarizes all results, among which the dimensionless residues, proposed as 
a guideline to choose the relevant mechanism complexity. Overall we find that 2-step series reaction 
models provide acceptable accuracy for reproducing the mass loss from lignocellulosic biomass and their 
constituents under torrefaction conditions. 

© 2013 Elsevier B.V. All rights reserved. 


1. Introduction 

Lignocellulosic biomass is considered as a potential renewable 
carbon neutral source of liquid fossil fuels. Various transformation 
processes exist, among which the BtL (biomass to liquid) process in 
which the first step requires the biomass to be pulverized into an 
EFR (entrained flow reactor) for gasification. The main difficulty is 
that lignocellulosic biomass is a fibrous and resilient material which 
is costly to grind and complicated to inject, two negative aspects 
regarding the efficiency of the BtL chain. 

Torrefaction is a mild thermal treatment at temperatures ran¬ 
ging from 200 to 300 °C under an inert atmosphere and with 
residence times of about 30-60 min. The process produces a mate¬ 
rial that has several advantages: it is easily ground to powders with 
good fluidization properties [ 1 ], it is hydrophobic and thus less sen¬ 
sitive to humidity upon storage [2,3] and it has increased LHV/HHV 
values [4,5], The main drawback is the overall energy loss of the 
process. A compromise must therefore be found between improved 
physical properties of the treated material and acceptable energy 
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loss. The extent of mass loss during torrefaction can be used as a 
first indicator of the alterations of biomass properties [ 5 ]. Since this 
loss is a strong function of the temperature history of the biomass, 
accurate kinetic models are a requirement for optimizing industrial 
operating conditions. 

Many kinetic models for biomass thermal treatment and for its 
main components (cellulose, hemicellulose and lignin) have been 
proposed in the literature with various degrees of complexity and 
are briefly reviewed hereafter. 


1.1. Models proposed for wood 

Early models [6] assumed the thermal treatment of wood to 
proceed via a single first order reaction giving a solid product and 
a gas phase. The mass loss rate is given by Eq. (1) where k(T) is 
the kinetic constant and follows the Arrhenius law as given in Eq. 
(2). The extent of conversion, a, is defined in Eq. (3) where X i0 is 
the initial mass of woody biomass, X, is the mass of the wood as a 
function of time and X in f is the mass of wood at the end of the exper¬ 
iment (thermally treated wood). In the case of isothermal tests, the 
hypothesis of a first order reaction (n = 1) allows Eq. (1) to be solved 
analytically to determine the rate parameters A, E and n, while in 
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tified using Eq. (4) I 


1 [?]■ 




Heating rate (K/min) Temperature (°C) E(kJ/mol) A(s-') 


10 

15 


40 


261-382 

260-402 

270-394 

280-406 

284-408 


117 

120 

122 

124 


3.32 x 10 9 
2.98 x 10 9 
6.68 x 10 9 
1.14 xlO 10 
1.9 xlO 10 


270-408 

272-425 

270-450 

278-458 

278-456 


120 


103 

103 


3.56 xlO 9 
1.26 xlO 9 
9.02 xlO 7 
1.51 x 10 s 
2.32 x 10 s 


the case of a dynamic tests with temperature ramps, the plot of the 
left hand term of Eq. (4) as a function of 1/T can be used. 

^=k(Tlf(«) = fe(T)(l-a)" (1) 


k(T)=Aexp^'j 


(2) 


j^iC 

X i0 


o-Xj 

.-Xinf 


(3) 


generates a solid and a gas phase. Later, Shafizadeh and Bradbury 
[13] modified this kinetic model with the hypothesis that the cel¬ 
lulose is first depolymerised or transformed into a more reactive 
phase (active cellulose) which is then transformed into three com¬ 
ponents: gas, tar and a solid residue. Currently, the most widely 
used model is the one proposed by Varhegyi et al. [14] (Table 2, 
scheme A) in which cellulose decomposes directly into gas and solid 
phases through two parallel reactions. Compared to the single first 
order reaction mentioned above (Eq. (1)) this kinetic representa¬ 
tion allows the variations of product yields with temperature to be 
taken into account. 


. (da/dt\ InA - E 

(4) 

Table 1 summarizes rate constant parameters identified by this 
method for thermal treatment of beech and pine [7] up to 400 °C 
with a constant heating rate varying from 10 to 40 K/min. The val¬ 
ues of E and A for beech are observed to increase with the heating 
rate, whereas the contrary is observed for pine. The variation of 
kinetic parameters with the heating rate when using a single step 
first order reaction has also been observed in coal pyrolysis [8] and 
is explained by the fact that coal devolatilization at the molecular 
level proceeds through a complex set of parallel and series reac¬ 
tions. Lignocellulosic biomass is expected to behave somewhat like 
coal in the sense that it also should evolve through many differ¬ 
ent reactions due to its complex molecular structure [9], Clearly a 
single-step reaction cannot be appropriate for describing the kinet¬ 
ics under all conditions. 

More recently Branca and di Blasi [10] proposed a three step 
kinetic model for wood thermal treatment in the temperature range 
of 265-435 °C as shown in Table 2 (scheme C). In this scheme, X 
is the wood, S is the thermally treated wood, V 2 and V 3 are 
volatiles, X' and X" are intermediate components. The transforma¬ 
tion was assumed to occur through three successive single first 
order reactions. Kinetic parameters were identified with several 
isothermal tests lasting 90 min after heating at 1000 K/min. 

With a similar approach, Prins et al. [ 11 ] modelled the mass loss 
of willow as a two step kinetic model under torrefaction conditions. 
As with Branca and di Blasi [10], all reactions were first order and 
followed the Arrhenius law. The kinetic parameters were identified 
under several isothermal tests at different temperatures, ranging 
from 230 to 300 °C with heat-up from 200°C at 10-100 K/min and 
observation times of 3 h. The two step mechanism was shown to 
provide good results above 270 °C but no comparisons were shown 
for lower temperatures. 

1.2. Models proposed for cellulose 

The first kinetic model for cellulose was developed by Broido 
and Weinstein [12] for temperatures higher than 400°C (Fig. 1). 
In the first step of this model, the cellulose generates some tar 
and anhydrocellulose, and in the second step the anhydrocellulose 


1.3. Models proposed for hemicelluloses 

Hardwood hemicelluloses are mainly composed of xylan (—80%) 
whereas softwood hemicelluloses are mainly made of glucoman- 
nan (—60%) [4], The latter is not as readily available commercially 
as the former, therefore studies on the thermal treatment of gluco- 
mannan are scare, and the literature mostly offers kinetic models 
for xylan. Just recently, Branca et al. [15] used a one and a two 
stage mechanism to describe the pyrolysis kinetics of commercial 
glucomannan. 

It has long been considered that xylan thermal treatment fol¬ 
lows an irreversible single-step first order reaction [16], However, 
as shown above for wood, such kinetic schemes lead to rate 
coefficients that are temperature and heating rate dependent. 
Multi-step series kinetic models are therefore more adapted for 
xylan weight loss prediction. Several two or three step kinetic 
mechanisms have been proposed [17-21] which provide more 
accurate predictions of hemicelluloses mass loss during pyrolysis. 

1.4. Models proposed for lignin 

To our knowledge, the only kinetic model existing for lignin tor- 
refaction is an irreversible single-step first order reaction. Several 
authors have identified kinetic parameters that vary according to 
the experimental conditions and the nature of lignin [17,22-25]. 

1.5. Summary 

In light of the above review, the various kinetic models pro¬ 
posed for predicting the mass loss of lignocellulosic biomass, xylan, 


Cellulose 


k. 



Anhydrocellulose 


Tar 


^3 0.65 Volatiles 

0.35 Char 


Fig. 1. Kinetic model of cellulose thermal treatment proposed by Broido and Wein¬ 
stein [12], 
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Table 2 

Summary of the main kinetic model existing on literature. 


Varhegyi et al. [14] 
Koufopanos et al. [17] 
Tang et al. [22] 
William et al. [23] 
Gronli et al. [24] 

Orfao et al. [25] 


Broido and Weinstein [12] 
Shafizadeh et al. [13] 
Varhegyi et al. [19] 
Lanzetta et al. [20] 

Fisher etal. [21] 


al. [10[ 
1. 





cellulose and lignin during thermal treatment, can be categorized 
into three types, A, B or C respectively for single, two or 3-step reac¬ 
tions as listed in Table 2. For example, the model proposed by Ward 
et al. [18] (Fig. 2) for the thermal treatment of xylan is considered 
to be one of type C. This model gives the proportion of volatiles for 
each irreversible first order reaction, with k v i /(k v i + k s i) = 0.16, 
W(kv2 + ks2) = 0.47 and k v3 /(k v3 + k s3 ) = 0.36. 

The kinetic parameters determined from these models have 
mostly been identified under dynamic tests with final temperatures 
above 400 °C exceeding typical torrefaction conditions. Thus it is 
not readily clear which type of model and complexity is required 
for adequately reproducing the rate of mass loss under moder¬ 
ate thermal treatment. Targeting conditions more representative 
of torrefaction conditions, new weight loss measurements are pre¬ 
sented here in the temperature range of 210-300 °C for various 
biomass and biomass component samples. The data are used to 
evaluate the different types of mechanistic schemes proposed in 


-► 0.16 M 1 


-► 0.84 Sj 




Fig. 2. Kinetic model for xylan thermal degradation developed by Ward et al. [ 18 J. 


the literature by using an inverse method to evaluate optimized 
kinetic parameters for each scheme and for each sample, valid 
for the whole range of temperatures studied. The most adequate 
mechanistic schemes are thus identified for each sample. 


2. Materials and methods 

2.1. Experiments 

Mass loss during torrefaction was measured by TGA under a 
nitrogen flow of 15 mL/min in a Setaram thermo-balance (Sensys- 
evo-3D) with less than 1% error. The accuracy of the measurements 
was estimated from several blank tests to correct the effects of 
buoyancy forces which depend on the gas flow rate and the tem¬ 
perature level. Each sample is first dried at 100° C for two hours 
then heated at 20K/min up to a final temperature of 210°C, 250°C 
or 290 °C which is maintained for a period of ten hours. This pro¬ 
tocol was designed to assess the kinetic rates under essentially 
isothermal conditions. 

Three species of wood have been studied: locust, a temperate 
hardwood, spruce, a softwood, and eucalyptus, a tropical hard¬ 
wood. They have been selected on the basis of their different 
structural properties which lead to different thermal resilience and 
thus possibly to different appropriate kinetic models. The weight 
loss of the main components of lignocellulosic biomass was also 
measured for cellulose (Sigma-Aldrich), xylan (Sigma-Aldrich), 
lignin (Green Value) and glucomannan (Kalys). The later is an 
important contribution of the present work as the literature is very 
poor in data regarding this important component of softwood. The 
polymers were supplied as powders, whereas the woody biomass 
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Fraction of mass loss that occurs during the heating phase. 


Temperature (°C) Fraction of mass loss at the end of the heating period (%) (AWL(t= lOh) (%)) 


Xylan Glucomannan Cellulose Lignin Spruce Locust Eucalyptus 


210°C 5.4 (26.1 ) a 0.8(11.2) «0(«0) 4(10.5) 1(8) 0.8(13.9) 1(9.3) 

250-C 11.6(51) 1 (51) 2.5(15.5) 11.8(17.2) 2.7(33) 2.7(33.3) 3.3(33.8) 

290-C 34.7(58.5) 3.6(66.2) 0.7(71.7) 13.6(24.8) 3.9(67.5) 6.1 (60.8) 6.5(58.9) 

a The first number is AW ^^ eall ' ls ■ 100, the number in parentheses is AWLioh (%). 


samples were ground to obtain particle sizes under 0.1 mm. Work¬ 
ing with powders allowed heat transfer limitations to be neglected. 
The sample initial weight varied from 15 to 18 mg for polymers and 
from 6 to 9 mg for woody biomass. 

2.2. Solution of the kinetics formulation 

For accuracy purposes, all kinetic models were solved numer¬ 
ically in this work. This allows any time-temperature pathway to 
be solved with the actual conditions undergone by the sample. The 
method principle will be detailed for model C. The strategy is simi¬ 
lar, and simpler, for models A and B (four differential equations for 
model B and only two for model A). 

In the case of model C described in Table 2, the set of differential 
equations describing the evolution of the system are written as 
follows: 


^r = -(ki+*2)|X] 

(7) 


(8) 

= fc 2 [X] - (fc 3 + k 4 )[X, ] 

(9) 

d[V 2 ] Ir ry 1 

T = k3[Xl] 

(10) 

d ^=k 4 lX x ]-(k 5 + k 6 )[X 2 ] 

(11) 


d[V 3 ] 

-dT = k5[X2] 

(12) 


(13) 

m=[X] + [X 1 ] + [X 2 ] + [S] 

(14) 


where [X] is the mass of the raw material, [Vi], [V 2 ] and [V 3 ] are 
the mass of volatiles, [S] is the mass of torrefied material, and [X] ] 
and [X 2 ] are the mass of intermediary components. All reactions are 
first order and kinetic rates k\, k 2 , k 3 and k 4 follow the Arrhenius 
law (Eq. (2)). 

If fully isothermal reaction conditions were reached the analyti¬ 
cal resolution of equations Eqs. (7)—(14) would have been possible. 
However, it is important to point out that the mass loss during the 
short heat-up phase is not systematically negligible. The ratio of 
this early mass loss relative to the total one observed after 10 h is 
specified in Table 3 for the various samples. The values are highest 
for xylan and lignin, exceeding 10% at temperatures above 250 °C. 

Due to this early mass release, the kinetic parameters need to use 
the actual sample temperature. Therefore, the set of Eqs. (7)—(14) 
was solved numerically, using a Crank-Nicholson scheme, a second 
order finite difference method in time, as described hereafter: 

d[X] ^ [X] n+1 - [X] n 
dt ~ dt 

& —0.5(k" +1 + k" +1 )[X] n+1 - 0.5(k" + k^)[X] n (15) 

” [Vlin+ dt~ 1Ul1 " ” 0.5fc? +1 [X] n+1 + 0.5k?[X] n (16) 



Fig. 3. Inverse method principle. 
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f ] Ri aC^w!,o.5( k f[xr 1 

-(fcf 1 + fcf 1 m ] n+1 ) + 0.5(k"[X] n - (k? + k n 4 )[X\ n (17) 

^ |%l n+ ^- [fig w 0 5(k n +1 [Xi jn+1 } + o.5(^[ Xl ]") (18) 

d\^] ^ [X2l " +1 d ~ [X2]n « 0.5(k^ +1 [Xi ] n+1 - (k" +1 + kg +1 )[X 2 ] n+1 ) 
+ 0.5(fc2[X 1 ] n - (k" + kg)[X 2 ] n ) (19) 

d[ £ ] « « 0.5k" +1 [X 2 ] n+1 + 0.5k" [X 2 ] n (20) 

w [ s l" +1 ~ PI" « 0.5kg +1 (X 2 ] n+1 + 0.5kg[X 2 ] n (21) 


(a) 210 °C 




In Eqs. (15)-(21), subscript n stands for time t and subscripts n +1 
stands for tine t+dt. All kinetic parameters k, are computed using 
the actual sample temperature (either time t or time t+dt), as 
obtained during the thermobalance tests. 

2.3. Parameter evaluation by inverse method 

Kinetic parameters were identified using an inverse method, the 
principles of which are summarized in Fig. 3. The originality of our 
approach is to fit a single set of parameter values using the three 
temperature levels simultaneously. A first set rate parameters with 
representative values are estimated and used to compute mass loss 
solutions for a given sample at the three temperatures. This compu¬ 
tation uses the experimental temperature histories as input data. 
An objective function is defined as the error between experiments 
and data over the full experimental times at all temperatures and 
the optimization algorithm varies the parameter values to mini¬ 
mize this objective function. The objective function F used in the 
present work is defined as follows. 


(d) 210°C 





Fig. 4. Measured anhydrous weight loss (symbols) of: (a-c) wood samples (spruce, 1 
lignin) for different temperatures (black lines). 


ucalyptus) and (d-f)' 


: (cellulose, xylan, glucomannan; 
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F = « 2 10^210 + «250f250 + “290^290 


EflKob 


(5) 


(6) 


and the otj are weighting coefficients for the different temperature 
levels (E a r = !)• N T is the size of the data vector at tempera¬ 
ture T, is the initial mass, m [ observed is the measured mass loss 
and m^ calcu | ate d is the corresponding calculated mass loss. As the 
total mass loss depends strongly on the temperature level (Fig. 4), 
the weighting coefficients aj have to depend on the temperature 
level to minimize simultaneously the global error and the partic¬ 
ular errors. The choice of these weighting coefficients was made 
to reduce the error at 210 °C while maintaining a good fit at the 
higher temperatures. A detailed description of the method derived 
to optimize the choice of the weighting coefficients is proposed in 
Section 3.2. 


3. Results and discussion 

3. 1. Weight loss experiments 

The anhydrous weight loss (AWL, defined as the ratio of the 
instantaneous mass loss as a function of time and the initial mass 
after drying) of spruce, locust and eucalyptus are compared at three 
different isothermal temperatures of 210°C, 250°C and 290°C in 
Fig. 4a-c respectively. Similar results for the lignocellulosic biomass 
constituents, cellulose, xylan, glucomannan and lignin are shown 
in Fig. 4d-f. 

The results display interesting relative behaviours between 
spruce, locust, xylan and glucomannan. At 210°C locust shows a 
higher AWL than spruce in agreement with the observations of 
others [11,25], This behaviour is explained by the different nature 
of hemicellulose in these samples. Locust hemicellulose is mainly 
composed of xylan, whereas spruce hemicellulose is mainly com¬ 
posed of glucomannan. As shown in Fig. 4a, at 210 0 C, spruce is more 
thermally stable than locust. However, for long residence times, 
reaching 10 h, and with increasing temperatures, the AWL of spruce 
close to that of locust at 250 °C (Fig. 4b) eventually supersedes it at 
290 °C (Fig. 4c). These relative characteristics are also observed for 
xylan and glucomannan which confirms the strong link between 
the reactivity of the wood and its composition in hemicelluloses. It 
should also be noted that softwood species, such as spruce or pine, 
are generally considered to exhibit lower weight losses than decid¬ 
uous trees. However, few studies extend to residence times of 10 h 
and our results show that at temperatures above 250 °C this trend 
eventually reverses. 

Finally, we note that lignin is the constituent that shows the 
least AWL variation with temperature and does not exhibit more 
than 20% AWL even at the highest temperature. In contrast, cel¬ 
lulose exhibits the strongest AWL loss with temperature (up to 


70%). As expected from the literature [4] cellulose exhibits very 
low mass losses at moderate temperature levels (below 250 °C), 
but presents the highest mass loss value at 290°C after 10 h of heat 
treatment. Nevertheless, the time constant for cellulose decompo¬ 
sition remains the largest one. 


3.2. Choice of the weighting coefficients 


The choice of the weighting coefficients otj is discussed here¬ 
after. Optimized kinetic parameters are found by minimizing the 
residue F. In principle we could have allowed free variation of the 
weight coefficients. However, this leads to poor results either for 
the higher or for the lower temperatures. Indeed, it was observed 
that kinetic parameters identified using only the test at 210°C, fits 
correctly the anhydrous weight loss measured at 210°C (1owF 2 io), 
but underestimates the anhydrous weight loss at higher temper¬ 
atures (high F 2 5 q and F 2 g 0 ). On the contrary, kinetic parameters 
identified at high temperature (250 or 290 °C) fit correctly the anhy¬ 
drous weight loss at high temperatures (low F 250 and F 290 ) but 
overestimate the weight loss at low temperature (high F 2 i 0 ). A 
constraint needs to be imposed to get a reasonable F T for each tem¬ 
perature (lower than 5%) while still having a low F value. A new 
dimensionless number is then introduced as described in Eq. (22) 
where Z is the ratio between the weighting factor at low tempera¬ 
ture (210 °C), aiow. and the weighting factor at high temperatures 
(250 and 290 °C), a high . 


Fig. 5a and b depicts the variation of the overall residue value (F) 
and of the residue values for each temperature level ( F T ) as a func¬ 
tion ofZ for xylan (a) and locust (b) in the case of B-type scheme. It is 
observed that the lowest value of Fis obtained forZ = 0, which is nor¬ 
mal, because in this case the identification by inverse method does 
not take into account measurements at 210°C (a2io = 0), which 
allows the identification to focus on levels 250 °C and 290 °C. Con¬ 
sistently, the error found for level 210°C is huge (more than 16% 
for the two examples depicted here). On the contrary, for large val¬ 
ues of Z, the overall criterion is rather good, the residue for 210 °C 
excellent, but those for 250 °C and 290 °C unacceptable. 

It is also observed that from Z= 6 for xylan and Z = 6.5 for locust, 
values of F 2 io, F 2 so and F 2 go become constant, while the overall 
residue value F is still decreasing. This is explained by the fact that 
for high Z values, the minimization algorithm practically does not 
take into account anymore the high temperature measurements. 
Indeed, as Z increases F decreases to a limit value which corresponds 
to F 2 jo value when F 2 so = F 29 o = 0 (or Z-» oo). 

By observing the curves obtained for either locust or xylan, an 
optimum appears clearly in the region of Z=2, for which both the 
global error and the individual errors remain below 5%. Therefore, 
the kinetic parameters presented hereafter have been evaluated 
with a 2 io = 0.5 and a 2 so = a 29 o = 0.25. 
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□ Measurement 210°C 
— Model B 210°C 
O Measurement 250°C 
—Model B 250°C 
» Measurement 290°C 
—Model B 290°C 
0 Measurement 270°C 
—Model B 270°C 


scales for xylan using model type B at different 



□ Measurement 210°C 
—Model B 210°C 
O Measurement 250°C 
—Model B 250°C 
x Measurement 290°C 
—Model B 290°C 
0 Measurement 270”C 
-Model B 270°C 


type B at different temperatur 


; (solid lines). 


3.3. Assessment of the identified kinetic parameters 

Typical results for xylan and locust are presented in Figs. 6 and 7 
in which the measured and predicted mass losses with a type-B 
model are compared. These figures prove that the identified kinet¬ 
ics are able to predict with good accuracy the measured mass loss of 
the samples throughout the entire experimental time, particularly 
for locust during the first hour, which includes the non isothermal 
heating phase (Fig. 7b), as well as up to ten hours (Fig. 7a). 

In order to validate the model with data, the kinetic models were 
tested on experimental data not used in the identification process. 
Results for the anhydrous weight loss of xylan and locust at 270 °C 
has also been measured in order to verify the model with a set of 
data that was not used in the kinetic parameters identification, the 
heating rate between 100°C and 270°C being also 20K/min. The 
computed mass loss using the experimental temperature history, 
according to kinetic parameters values identified with measure¬ 
ments at 210,250 and 290 °C, is found to agree well with this data 
at 270 °C. 


Finally, as further test of the quality of the proposed kinetics 
model, Fig. 8a and b represents the comparison between the mea¬ 
sured and predicted mass loss rate (d(AWL)/dt) for locust and xylan, 
respectively, as a function of time with the type-B model. First, it 
is observed that a double peak of mass loss rate appears for xylan 
during thermal treatment at 290 °C, while no visible double peak 
exists during thermal treatment of locust. However, model B is able 
to predict with an acceptable accuracy the mass loss rate even in 
the case of two maximal values of mass loss rate. 

The identified value of kinetic parameters for all samples are 
presented in Tables 4-6 for types A, B and C models respectively. 
For the wood samples, in the case of type B and C models, the first 
reaction step is fast and is not rate limiting. This result is in agree¬ 
ment with the findings of Lanzetta et al. [20] and Prins et al. [11 ], 
and is coherent with the fact that the first step should represent 
hemicelluloses decomposition whereas the last one accounts for 
cellulose decomposition. 

The overall residue values obtained with models A, B and C are 
presented in Fig. 9 for all samples. It can be seen that models B and 















































■al/Th, 



Xylan 

Cellulose 

Glucomannan 

Lignin 

Locust 

Eucalyptus 

Spruce 

At (s~') 

El (kj/mol) 

A2 (s _1 ) 

E2 (kj/mol) 

2.9 xlO 9 
129.9 

2.9 x 10 5 
89.2 

6.3x10 17 

230.1 

4.2 x 10 11 

168.1 

3.9 xlO 10 

147.8 

4.1 x 10 7 

118.5 

2.4 xlO 7 

114.8 

3.1 x 10= 

88.7 

4.3 x 10 5 

102.2 

0.15 

34.9 

2.8 x 10= 

3.5 

48.9 

4.4 x10 s 

136.9 

1.78 x 10 s 

135.9 

r, s , 









Xylan 

Cellulose 

Glucomannan 

Lignin 

Locust 

Eucalyptus 

Spruce 

A1 (s" 1 ) 

El (kj/mol) 
A2(s-') 

E2 (kj/mol) 

A3 (s’ 1 ) 

E3 (kj/mol) 

A4 (s -1 ) 

E4 (kj/mol) 

2.2 xlO 9 
121.7 

2.7 xlO 2 

91.7 

66.6 

4.1 x 10 4 

89.5 

1.2 x10 s 

233.8 

6.5 xlO 9 3) 

157.3 

28.3 

33 xlO 3 

80.3 

1.7x10" 

<10 7 

114.9 

2.7x10= 

87 

1.9x10= 

99 

3.45x10= 

102.9 

1.8x10= 

81.4 

1.24 xlO 4 

96.9 

9.81 x 10 3 

85.9 

4.1 x 10 9 

133.9 

1.35 x 10 7 
101.7 

9.9x10= 

5.3 xlO" 3 

213 

9.3 x 10 7 

119.2 

2.4 xlO 3 

64.9 

4.3 x 10= 

106.4 

5.2 

52.3 

1.4 xlO 7 
113.2 

1.5 x 10 3 
62.9 

1.6x10= 

112.6 

7.7 x 10- 3 
35.5 

Table 6 









Xylan 

Cellulose 

Glucomannan 

Lignin 

Locust 

Eucalyptus 

Spruce 

A1 (s') 

El (kj/mol) 
A2(s-') 

E2 (kj/mol) 
A3(s-') 

E3 (kj/mol) 
A4(s-') 

E4 (kj/mol) 
A5(s-') 

E5 (kj/mol) 

ssu 

2.06 xlO 9 
120.4 

53 

43 

2.36 xlO 2 

70.3 

4.09 xlO 4 

88.3 

1.42 x 10 4 

738 x 10 2 

48.4 

1.29 xlO 19 
244.4 

6.93 xlO 9 

145.1 

34.02 

75.8 

8.44 xlO 3 

78.7 

7.9 xlO 4 

93.7 

5.79x10= 

101.2 

3.58 xlO 12 
166.3 

1.19 xlO 7 

108.5 

488.34 

74 

2.95x10= 

97.9 

1.75 xlO 4 

2X17 xlO 4 

2.6 xlO 7 

2.2 xlO 4 

76.3 

1.3 xlO 4 

86.2 

1.2 xlO 4 

75.1 

9.6 xlO 3 

97.7 

4.9 10 3 

84.5 

3.9 xlO 9 
131.9 

1.2 xlO 7 
100.2 

9.7 x 10= 

6.3 x 10- 3 
20.3 

1.3 xlO 4 

1.5 xlO 4 

90.1 

1.7 xlO 7 

112 

628.2 

58.7 

1.1x10= 

110.1 

15.9 

56.4 

1x10= 

4.2 x 10 4 
513.4 

1.2 xlO 7 
111.4 

1.6 xlO 3 
62.5 

1.5 x 10= 
113.1 

7.8 xlO" 3 
37.9 

5.1 x 10= 

6.3 x 10 4 

78 



Fig. 9. F values for kinetic model A, B and C for all samples with Z=2. 


C always give a lower residue than model A which means that mod¬ 
els B and C allow a better prediction of the mass loss than model 
A. In addition, model C does not significantly decrease the residue 
compared to model B, except for lignin. Therefore model B is rec¬ 
ommended for xylan, glucomannan, cellulose, spruce, locust and 
eucalyptus and model C for lignin. 

In the case of cellulose, one can see that the residue with model 
types B and C remain high compared to those obtained for the 
others samples. The main difficulty is that these kinetic schemes 
do not allow a proper description of cellulose depolymerisation 
when the treatment temperature approaches 300 °C. Some form of 
parallel reaction step is possibly required, and more sophisticated 
mechanism should be developed to improve predictions. 


4. Conclusion 

The anhydrous weight loss (AWL) during isothermal treatment 
under inert atmosphere of three woody biomasses (locust, spruce 
and eucalyptus) and their main components (xylan, glucomannan, 
lignin and cellulose) have been measured using an accurate ther¬ 
mobalance (Sensys-evo 3D). We observed that the higher mass 
loss rate observed in locust compared to spruce can be partially 
explained by the higher reactivity of xylan (main component of 
locust hemicelluloses) relative to that of glucomannan (main com¬ 
ponent of spruce hemicelluloses) only for moderate treatment 
intensity. At high temperatures (290 °C) and long residence times 
(5 h) spruce decomposition is stronger due to both higher cellulose 
and glucomannan reactivities. 

Based on this comprehensive set of experimental data, kinetic 
models suitable for implementation into a macroscopic computa¬ 
tional model [26,27] have been evaluated. The kinetic parameters 
for multi-step series reaction mechanisms found in the literature 
(1,2 and 3 step-models) have been evaluated by an inverse method. 
This method computes the solution of the multi-step formulation 
using the experimental temperature histories undergone by the 
samples as input data. In addition, the inverse procedure takes into 
account simultaneously the three temperature levels for all sam¬ 
ples, resulting in a prediction model valid in the whole studied 
temperature range. 

The results show that the accuracy of a 2-step series reaction 
model is generally satisfactory for representing the mass loss of 
lignocellulosic biomass and their components. The gains from a 
more complex 3-step model are not sufficient to warrant the added 
complexities. It therefore appears that further improvements in 
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thermal decomposition models in this temperature range would 
require other kind of kinetic schemes. 

Further work is currently undertaken to test the potential 
of these formulations when implemented in our comprehensive 
heat and mass transfer computational model [28], both at the 
micro-particle scale at different temperatures as well as at the 
macroscopic scale of a fixed bed reactor. 
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